Gravitational radiation from the r-mode instability 



Benjamin J. Owenf, and Lee LindblomJ 

f Department of Physics, University of Wisconsin-Milwaukee, P.O. Box 413, 
Milwaukee, WI 53201 USA 

I Theoretical Astrophysics 130-33, California Institute of Technology, 
Pasadena, CA 91125 USA 

E-mail: owen@gravity . phys . uwm . edu, lindblom@tapir . caltech . edu 

Abstract. The instability in the r-modes of rotating neutron stars can (in principle) 
emit substantial amounts of gravitational radiation (GR) which might be detectable 
by LIGO and similar detectors. Estimates are given here of the detectability of this 
GR based the non-linear simulations of the r-mode instability by Lindblom, Tohline 
and Vallisncri. The burst of GR produced by the instability in the rapidly rotating 
1.4 Mq neutron star in this simulation is fairly monochromatic with frequency near 
960 Hz and duration about 100 s. A simple analytical expression is derived here for 
the optimal S/N for detecting the GR from this type of source. For an object located 
at a distance of 20Mpc we estimate the optimal S/N to be in the range 1.2 ~ 12.0 
depending on the LIGO II configuration. 
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1. Introduction 

This paper contains a portion of the material presented by Lee Lindblom in the talk 
"Relativistic Instabilities in Compact Stars" at the Fourth Amaldi Conference on 
Gravitational Waves. Most of the material presented in that talk has already been 
published, and so we limit our discussion here to that material which has not previously 
appeared in print. We direct those readers who might be interested in the wider range 
of subjects covered in the talk to the recent review Lindblom (2001), and the recent 
papers Lindblom, Tohline and Vallisneri (2001b) and Lindblom and Owen (2001). 

Gravitational radiation (GR) is a de-stabilizing force on the r-modes of rotating 
neutron stars (Andersson 1998, Friedman and Morsink 1998). The first estimates of 
the timescales associated with this instability (Lindblom, Owen and Morsink 1998) 
showed that the GR driving force is much stronger than the stabilizing effects of the 
simplest forms of viscous dissipation in neutron stars. As a consequence there has 
been a great deal of interest in this instability as a potential source of observable GR, 
and as a mechanism for removing angular momentum from rapidly rotating neutron 
stars. During the past several years the effects of a number of additional dissipation 
mechanisms on the r-mode instability have been studied in some detail. A number of 
these mechanisms are much more effective in suppressing the r-mode instability than 
the simple viscosity considered in the initial estimates. In particular the effects of a 
solid crust (Bildsten and Ushomirsky 2000, Lindblom, Owen and Ushomirsky 2000, 
Wu, Matzner and Arras 2001), the effects of magnetic fields (Rezzolla, et al 2001a, 
2001b, Mendell 2001), the non-linear effects of mode-mode coupling (Schenk, et al 2001), 
and the effects of hyperon bulk viscosity (Jones 2001a, 2001b, Lindblom and Owen 
2001) make it appear less likely that the GR instability in the r-modes will play an 
interesting role in astrophysics. However, at the present time none of these mechanisms 
is understood well enough for us to conclude absolutely that the r-mode instability will 
never play a role in any neutron stars. Thus for the purposes of the present paper, 
we assume that the instability will occur in some rapidly rotating neutron stars. Our 
aim here is to present the best estimate of the GR that might be emitted during such 
an instability. To do this we analyze the GR emitted by the best currently available 
numerical simulation of the non-linear evolution of the r-mode instability. This paper 
is an update of the initial estimates of the GR emitted by the r-mode instability given 
by Owen, et al (1998). 

2. r-Mode Evolution Model 

We base our estimates of the GR produced by the non-linear evolution of an unstable r- 
mode, on the numerical simulations by Lindblom, Tohline and Vallisneri (2001a, 2001b). 
The evolution of a small amplitude m = 2 r-mode subject to the current-quadrupole 
GR reaction force was studied numerically in a 1.4M polytropic neutron star rotating 
at about 95% of its breakup angular velocity. Figure [l| illustrates the evolution of the 
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Figure 1. Evolution of the current quadrupole moment J22 ( m cgs units) from the 
numerical simulation of the GR driven growth of the m = 2 r-modc. 



current quadrupole moment | J 22 1 (in cgs units) of this model as the simulation evolves. 
The current quadrupole moment is defined as 



type vector spherical harmonic (Thorne 1980). The time displayed along the horizontal 
axis of Fig. p] is given in units of the initial rotation period of the star: Pq = 1.18 ms. 
During the first part of the evolution GR drives the growth of the r-mode, leading to 
the exponential growth in J 22 illustrated here. Once the amplitude of the mode becomes 
sufficiently large, however, non-linear hydrodynamic processes also become important. 
These lead to large surface waves which break and shock (see Lindblom, Tohline and 
Vallisneri 2001a, 2001b). The dissipation in these shocks damps the r-mode and leads 
to the rapid decrease in J 22 following its peak. 

The timescale of the GR instability in the r-modes is much much longer than the 
hydrodynamic timescale. To perform the numerical simulation illustrated in Fig. [1] it 
was necessary to increase artificially the strength of the GR driving force so that the 
instability could proceed more rapidly. Fortunately tests have shown (see Lindblom, 
Tohline and Vallisneri 2001b) that the maximum amplitude which the r-mode achieves 
(and hence the maximum value of J 22 ) is relatively insensitive to the strength of the GR 
driving force. The amplitude grows until the velocities in the mode reach some critical 
value before the waves break and shock. The value of this critical velocity appears to be 
relatively insensitive to how rapidly the fluid velocity is increased to this critical level. 

To determine the GR signal that would be emitted by the growth of an unstable 
r-mode, it is necessary to re-scale the time in the simulation. The early part of the 
simulation has been designed to proceed more rapidly than the physical case by the 
factor TG#(physical)/rGi?(simulation), where tqr represents the GR growth timescale of 
the instability. This ratio has the value 4488 in the simulation used here. Thus the early 
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Figure 2. Evolution of the "frequency" of the r-mode. The time parameter used 
here (in seconds) is scaled from the simulation to reflect the physical case until the 
formation of shocks. Noise early in the simulation is due to contamination of the 
current quadrupole J22 by other modes. 



exponential growth phase of a physical r-mode evolution will last longer than this phase 
of the simulation by this factor. The late stages of the simulation are dominated by 
hydrodynamic forces which damp the r-mode within a few rotation periods. This final 
stage of the evolution therefore proceeds at the same rate in both the physical and the 
simulation cases. We (somewhat artificially) choose the dividing time between the early 
and late stages of the evolution to be the time at which the value of | J22I is maximum. 

The time dependence of J22 is found in the simulation to be nearly sinusoidal with 
a time dependent amplitude, J22 = |^22| e ^ > (with dtp/dt = u ~ constant) because 
J22 is dominated by the unstable r-mode in this case. The time dependence of | J22I is 
illustrated in Fig. |l|. The frequency of the sinusoidal time dependence is conveniently 
determined numerically using the formula 



This approximation introduces errors of order {ujtgr}^ 1 ~ 2% in our simulation. 
Figure |2| illustrates the evolution of the frequency, / = — 2ttlu, determined numerically 
in this way. Here we plot the frequency as a function of the re-scaled physical time (in 
seconds). We note that the frequency changes by only a few percent during the course 
of the evolution in which about 40% of the angular momentum of the star is radiated 
away as GR. 

In general relativity theory an isolated object with time dependent current 
quadrupole moment radiates GR. The expression for the dimensionless GR amplitude 
h (averaged over possible source and detector orientations) is given by 




(3) 
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Figure 3. Evolution of the angle averaged dimensionless gravitational wave amplitude 
\h(t)\ for a source located at a distance D = 20Mpc. The time parameter used here (in 
seconds) is scaled from the simulation to reflect the physical case until the formation 
of shocks. 



where G and c are Newton's constant and the speed of light, and D is the distance to the 
source. Figure [5] illustrates the time dependence of this dimensionless GR amplitude as 
a function of time as determined by the simulation. Since the frequency of these r-mode 
oscillations is essentially constant, the GR amplitude just grows in proportion to J 22 . 
Once the amplitude peaks, non-linear hydrodynamic forces (i.e. shocks) quickly damp 
the mode and this is reflected in the sharp drop in h (as a function of physical time) 
illustrated in Fig. |[ We have chosen D = 20 Mpc, the distance to the Virgo Cluster 
of galaxies, for the purposes of the illustration. Optimistic estimates of the event rate 
for these objects suggests that the nearest events expected during the time frame of the 
LIGO observations will be at this distance. 



3. Optimal S/N 

The optimal signal to noise ratio (S/N) for detecting a GR signal can be achieved with 
the use of an optimal filter, consisting of a template that matches the waveform of the 
signal. Let h denote the Fourier transform of the GR signal: 

_ POO 

h(f) = / h(t)e- 2nift dt. (4) 



The discrete Fourier transform of the gravitational wave signal from our numerical 
simulation is shown in Fig. f|. Here we have smoothed the resulting \h(f)\ 2 with a 
windowing function of width 0.5 Hz. We see that the spectrum of the gravitational 
radiation is confined essentially to the band 940 < / < 980 Hz, which corresponds well 
with the evolution of the frequency of the r-mode as shown in Fig. |^. For a complex 
signal h(t) the optimal value of S/N is given by 

s\ 2 9 r \Kn\ 2 df ,,, 
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Table 1. Coefficients for the Taylor expansion of the one-sided power spectral density 
Sh at / = 960 Hz for three plausible LIGO II configurations. 
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Figure 4. Fourier transform of the gravitational wave amplitude illustrated in Fig. H, 

~ 1/2 

\h(f)\. Also shown are the noise power spectral density S h for the three potential 
LIGO II configurations used in our estimates of the optimal S/N. 



where Sh{f) is the one-sided power spectral density of detector noise. (For a real signal 
the 2 becomes a 4.) We approximate Sh{f) for frequencies near 960 Hz by the Taylor 
expansion 

S h (f)^S + S'Af + 0.5S"(Af) 2 , (6) 

where Af = f — 960 Hz, and the three constants S, S' and 5"' are listed in Table |I] 
for three plausible LIGO II design options (see Buonanno and Chen 2001). Using these 
values for Sh we find S/N w 1.2, 4.0 and 10.4 by performing the integral in Eq. (|5|) 
numerically for these three possible LIGO II configurations. 

We can also make an analytical estimate of the optimal S/N using an extension 
of a very general argument given originally by Blandford (1984, unpublished). For GR 
emitted by a multipole with azimuthal quantum number m the angular momentum loss 
rate is 

f = -^/^>W| 2 , (7) 

where |/i(t)| 2 has been averaged over all possible orientations of source and detector 
(and over several wave periods for the case of a real signal). When a function such as 
h(t) involves a rapid oscillation together with a much slower evolution of its amplitude 
and frequency, the Fourier transform h(f) is well approximated by the stationary phase 
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approximation, 



= \h{f)\ 2 | (8) 

for a complex signal h(t). (For a real signal, the left-hand side is averaged over several 
periods and the right-hand side is multiplied by 2.) Since the angle-averaging affects 
both sides of Eq. (H) equally, we can use Eqs. (0) and @ to re-express S/N for a source 
at an "average" orientation from Eq. (|5]) as 

S\ 2 AG r dJ 

(9) 



MJ 5m7rc 3 D 2 J fS h (f)' 
For GR sources such as the r-mode evolution the frequency of the radiation emitted is 
nearly constant, so fSh(f) can be treated as being essentially constant in Eq. (||). Thus, 
this integral becomes just the total amount of angular momentum |AJ| radiated away 
as GR: 

S\ 2 _ AG \AJ\ 

.NJ ~ hrrnic^D 2 fS h {fy ( ' 

The total amount of angular momentum radiated away as GR in the simulation was 
|AJ| 4.5 x 10 48 in cgs units. Using this value in Eq. (|10D together with the values of 
Sh from Eq. (||) and Table [I], we find S/N « 1.4, 4.5, and 12.0 in good agreement with 
our previous estimate based on the discrete Fourier transform and direct integration of 
the simulation waveform. We note that this argument can easily be extended to signals 
composed of multiple harmonics and multipoles so long as their frequencies are well- 
defined: Time- averaging of Eq. (^) eliminates cross-terms, so Eq. ( ]TD| ) is simply summed 
over each harmonic and m. We can also extend the argument to non-monotonically 
evolving frequencies by using Eq. (^j) to express Eq. (0) as an integral over time and 
dividing f(t) into piecewise monotonic parts. 

It is not unreasonable to think that a realistic data analysis strategy could come 
within a factor of 2 of the optimal S/N. Due to the complexity of the physics involved 
it seems unlikely that matched filtering will ever be a viable option. However, cross- 
correlation of the output of two aligned interferometers (LIGO Hanford and LIGO 
Livingston) can in principle achieve l/v2 of optimal S/N if the detectors' noise is of 
comparable strength and uncorrelated — including non- Gaussian noise bursts (Anderson 
et al 2001). This strategy relies on the supernova associated with the r-mode having 
been observed optically, allowing the appropriate time delay to be inserted between 
the interferometer data streams. Narrowing the search to a few minutes after the 
supernova (instead of continuous operation) also has the effect of greatly reducing 
the S/N threshold for detection with reasonable false alarm statistics. Thus a cross- 
correlation with S/iV ~ 4 might be considered enough for detection, implying a 
realistically detectable distance for r-modes of 5 Mpc for even the least optimal (for 
this type of source) LIGO II configuration — and up to 50 Mpc for the most optimal. 
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